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Abstract 

Recently we studied the effects of information carrying waves prop- 
agating through isotropic cosmologies. By information carrying we 
mean that the waves have an arbitrary dependence on a function. We 
found that the waves introduce shear and anisotropic stress into the 
universe. We then constructed explicit examples of pure gravity wave 
perturbations for which the presence of this anisotropic stress is essen- 
tial and the null hypersurfaces playing the role of the histories of the 
wave-fronts in the background space-time are shear-free. Motivated 
by this result we now prove that these two properties are true for all 
information carrying waves in isotropic cosmologies. 
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1 Introduction 



In a recent paper £Q the propagation of shear-free gravitational waves 
through isotropic universes was studied using the gauge-invariant and covari- 
ant formalism of Ellis-Bruni [2|. In this approach the waves are represented 
as a perturbation of the Robertson-Walker space-time. The gravitational 
waves studied also had the property that they could carry arbitrary informa- 
tion. By this we mean that the Ellis-Bruni gauge-invariant variables had an 
arbitrary dependence on a function. 

It was found that the consistency of the partial differential equations gov- 
erning the perturbations requires that all gauge-invariant perturbed quan- 
tities must vanish except for the perturbed shear of the matter world-lines 
and the anisotropic stress of the matter distribution. A wave equation for 
the perturbed shear was derived. In addition the perturbed shear was shown 
to satisfy a propagation equation along the null rays associated with the 
gravitational waves (see eqs. (|2.7J1 and ()2.8|) below). 

For pure gravitational radiation (i.e. the perturbed Weyl tensor is type N 
in the Petrov classification) explicit examples were constructed for which it 
was shown (using the wave equation and the propagation equation) that the 
presence of the anisotropic stress is essential and the hypersurfaces playing 
the role of the histories of the wave-fronts are shear-free. We now prove that 
these two properties are true in general and are not in any way dependent 
on the examples constructed. 

The paper is organized as follows: The notation used and some useful 
equations are given in the Appendix and referred to where appropriate. In 
Section 2 we outline the results obtained in |X| which we use in the paper. The 
first of our results is also given in this section. The elements of the Newman- 
Penrose formalism necessary for this problem are described in Section 3 and 
the main result is derived in Section 4. The paper ends with a summary of 
the results of Sections 3 and 4. 

2 Equations Governing Gravitational Waves 

Throughout this paper we shall use the notation and sign conventions of 
[H] . For convenience these have been briefly outlined in the Appendix. We are 
concerned with a four dimensional space-time manifold with metric tensor 
gij in a local coordinate system x l and a preferred congruence of world-lines 
tangent to the unit vector field u l (with u % Ui = —1). The decompositions 
of the symmetric energy-momentum-stress tensor T^- and the Weyl tensor 
with respect to this 4-velocity field, along with definitions of the kinematical 
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quantities cry, CJy, 6 1 , can be found in the Appendix. 

In a recent paper [I] we used the gauge-invariant and covariant ap- 
proach of Ellis-Bruni 2J to construct gravitational wave perturbations of 
Robertson-Walker space-times. The background metric tensor is thus the 
Robertson- Walker metric, the background energy-momentum-stress tensor 
is (jA.9|) specialized to a perfect fluid (i.e. with q l = = tt 1 - 7 ) with fluid 
4-velocity u l and the background Weyl tensor vanishes. The Ellis-Bruni 
approach involves working in a general local coordinate system with gauge- 
invariant small quantities rather than small perturbations of the background 
metric. These quantities have the property that they vanish in the back- 
ground space-time. For isotropic space-times the Ellis-Bruni variables are 
dij, u\ uJij, Xi = hlfij, Y { = hfpj, Zi = hjOj, Try, ft, £y and #y (see 
the Appendix for definitions of these variables). We found that it is tensor 
quantities that describe gravitational wave perturbations and thus we can 
set all Ellis-Bruni variables equal to zero except for Try, oy, £7y and f/y. 
The equations satisfied by these variables are obtained by projections, in the 
direction of and orthogonal to the 4-velocity u l , of the equations of motion 
and the energy conservation equation contained in T^.j = 0, the Ricci iden- 
tities and the Bianchi identities (see Appendix). These projections give rise 
to a lengthy list of equations which can be found in pQ. 

We specialized to a particular class of gravitational waves, namely those 
which can be viewed as carrying arbitrary information. To incorporate this 
into the Ellis-Bruni set up we required that the gauge-invariant small quan- 
tities have an arbitrary dependence on a function. Specifically we assume 
that 

Try = Ily F{<t>) , on = syF(0) , (2.1) 

where F is an arbitrary real-valued function of its argument 4>{x l ). This 
idea of introducing arbitrary functions into solutions of Einstein's equations 
describing gravitational waves goes back to work by Trautman |3] and this 
form for the Ellis-Bruni variables was first used by Hogan and Ellis [H] . Using 
the projections of the Ricci identities with uj^ = = u l we can write 

12 2 
E v = 2 n v + 3 °" 2/ ^ ~ 3 Ui i ~ ° ih ak 'j ~ h i h 3 &kl ' ( 2 - 2 ) 

and 

H ij = -h*h l j a ih ™r h )r Pq u r . (2.3) 

Thus these variables are derived from and TTy and it is not necessary to 
make any extra assumptions about their dependence on F(4>). We note that 
Hij, Sij are orthogonal to u % and trace-free with respect to the background 
metric ftj. 
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Substituting (j2.1|) into the linearized versions of the equations satisfied 
by these variables we find that and IT J are orthogonal to the gradient of 
and divergent free i.e. we have 

«' J o.,, = , II'-'o.; = , (2.4) 

and 

8*3,3 = , U'r , = . (2.5) 
Also, with <f> := 4>jU l ^ 0, consistency of the equations requires that 

ll'o^o.j 0, (2.6) 

where is the background metric. This means that the hypersurfaces 
(j)(x l ) = const, in the background isotropic space-time playing the role of 
the histories of the wave-fronts must be null. The following wave equation 
for s y ', 

s ir ' k -,k + $ ° 2 ) slJ + (p-\ f°) s ij = -n iJ - \ e ir j ' , (2.7) 

and propagation equation for s i3 - along the null geodesies tangent to 4>'\ 

Sir,k<P' k + (g^*:*- 3 d< j>) s ij = ~2^ii » ( 2 - 8 ) 

are consequences of the Ricci identities. We have used fio here to denote the 
energy density of the cosmic fluid to avoid confusion with the spin coefficient 
/i used in the next section. The internal consistencies of these equations were 
checked in pQ. Substituting (|2.1j) into ()2.2j) and ()2.3)1 we find that the electric 
and magnetic parts of the Weyl tensor are now given by 

Eij = (- U i3 - - | 9 s tj ^j F-4> Sij F' , (2.9) 

and 

H tj = -s ( /' m Vj )ki m u k F - s { / r] j)klm u k 0' m F' , (2.10) 

where F' — 

We are interested here in pure gravitational wave perturbations i.e. per- 
turbations which have pure type iV perturbed Weyl tensor in the Petrov 
classification. It is easily checked using Eq. (|231) that the F -parts of E {j 
and Hij above are automatically type N with degenerate principal null di- 
rection c()' 1 . However the F-parts of and are not in general type N 
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and so to describe pure gravitational wave perturbations we require (see [I]) 
that 

s i3 >,i = 0, - ! o., :i '-/o, : ' ! . (2.11) 

Making use of the following null tetrad, ki = — (p^ 1 0^, k — U{ — hki, and 
rrii, fhi a complex covariant vector field and its complex conjugate chosen 
so that they are null (m 1 = = m l rhi), are orthogonal to k 1 and l l and 
satisfy rn'm; = 1 we find that (see [T]) the conditions (|2.11|) are equivalent 
to 

s^i-jfrfP = (2.12) 

and 

s (f) ^. jim 1 m 3 = s (f) t i-jfh l fnP . (2-13) 

In terms of the tetrad we can write s u as (because s lJ Uj = 0, s l j = 0, 
s lj (f)j = and so s l %j = and s l Hj = 0) 

s'J = sm l m j + sfh l fh j , (2.14) 

where |s| 2 = xs^'sjj. For future use we point out here that IP 7 can be 
expressed on the tetrad in a similar way, namely, 

U ij = 5mW + nmW . (2.15) 

We emphasize that we are using II here for the complex component of Il afe 
to avoid confusion with the spin coefficient tt appearing in the next section. 
Taking the divergence of the first of Eqs. (j2.4j) and substituting for s iJ from 
tiTTty yields (since s ij .j = 0) 

(srrfm? + sm'm^^j.j = . (2-16) 

Comparing this with Eq. ()2.13|) it is clear that for consistency (with s ^ 0) 
we must have 

hi-ynfm? = . (2.17) 

Thus the hypersurfaces 4>(x l ) = const, which play the role of the histories of 
the wave-fronts must be shear-free. With s ^ ()2.12|) becomes 

0^.77^ = 0. (2.18) 

This condition corresponds to the vanishing of the spin coefficients r, v which 
we shall make use of in the next section. 
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3 The Newman— Penrose Formalism 

For the remainder of the paper we seek to demonstrate that the presence 
of anisotropic stress in the perturbed space-time is essential for the existence 
of pure gravitational waves. We find the most effective way to do this is to 
make use of the Newman-Penrose Formalism. A comprehensive description 
of this formalism can be found in [Bj. Specifically we choose the tetrad 
described after Eq. (|2.11j) and we label the vectors as follows: 

D = k* di = e\ di , A = P di = e\ d t , , , 

5 = m l di = e\di , 5* = fn l di = e\ di . ^ ' ' 

The spin coefficients (Ricci rotation coefficients) are denoted by 7 a b c = — 7b ac - 
It follows from Eq. (j2.17|) that k 1 is shear-free (i.e. ki-jirfmP = 0) and a 
simple calculation using the fact that we can write tp^ = (</>,&)' + ^9hl<j)^ c 
shows that k l is also geodesic with hf-jti = ^9k\ As a consequence of these 
properties the spin coefficients k, a, A and tc vanish. We again point out 
that 7r here is not to be confused with the complex component II of IF 3 . 
As previously mentioned after Eq. (|2.18J) the spin coefficients v and r also 
vanish. Therefore, for the problem at hand the only non-zero background 
spin coefficients are p, fx, 7, e, a, (3 with 



p = 




(3.2) 


p = 


w 

T 3 ' 


(3.3) 


e = 


-(7211+7341) , 


(3.4) 


7 = 


l , \ 
2(7212 + 7342) , 


(3.5) 


a = 


-(7214 + 7344) , 


(3.6) 


P = 


!/ \ 

2 (7213 + 7343) • 


(3.7) 



Using 7 a {, c = — 7b ac and noting that the 4-acceleration vanishes in the 
background we find that 

7 + 7* = ^, e + e* = ~9, (3.8) 
o 3 

where, following Chandrasekhar, the star denotes complex conjugation. In 

addition is a real-valued function and thus p and p are real and we have 
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p-p* = 0, fi-fi* = 0. (3.9) 

We will use these facts in the next section when we examine the wave equation 
()2.7|) and the propagation equation ()2.8j) in this formalism. 
In the tetrad formalism the Ricci identities ()A.7j) read 



Rabcd = -Jabc,d + labd,c + Jbafilc d ~ Id c) + IfacW d ~ IfadTfl > c , (3.10) 

where R a bcd are the tetrad components of the Riemann tensor. An alternative 
expression for these components is given via the definition of the tetrad 
components of the Weyl tensor, namely, 

Rabcd = C abcd + \ ( Vac R bd - Vbc Rad ~ VadRbc + VbdRac) 

~ ■ziVacVbd-VadVbcjR ■ (3.11) 


Here R ab = Rije a e{ are the tetrad components of the Ricci tensor, R = /x — 3p 
is the Ricci scalar and r\ ab := e\en, = except for r/i 2 = 7721 = — r/34 = —^43 = 
— 1. In general we label tetrad components of tensors using a, b, c, . . . and 
coordinate components using k, . . . .Making use of Eqs. (jA.8|) and (jA.9|) 
we find that the tetrad components of the background perfect fluid can be 
written in the form: 

Rab = (A*Q + P)U a U b +^(^0-p)Vab , (3. 12) 

where u a = Uit a with Ui = U + \k{ and thus u a = (— 1,— 1,0,0). Here 
u a are the tetrad components of the 4-velocity of a fluid particle, /io is the 
proper density and p is the isotropic pressure. We shall always assume that 
/i + p ^ so that R ab has a unique time-like eigenvector. Evaluating 
the various components of the Riemann tensor using both Eqs. (j3~Tnj) and 
()3.11|) (we use the commutation relations involving the operators D, A, 5, 5* 
to do this) and then equating the two different expressions obtained yields 
the following twelve equations for the derivatives of the background spin 
coefficients: 

Dp = -p 2 -(e + e*)-^o+p) , (3.13) 
Da - S*e = -a(p + e* - 2e) + (3*6 , (3.14) 
DP - 5e = p(e* - p*) + ea* , (3.15) 
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D 1 -Ae = 2e 1 + e 1 * + e* 1 --fi --p, (3.16) 

Dfi = -iip* + /x(e + e*) + -p- —fio , (3.17) 

Sp = -p( a * + (3) , (3.18) 

5a - 5* (3 = -up - aa* - (3(3* + 2a(3 - -/i , (3.19) 

6 

5*n = p(a + (3*) , (3.20) 

A/i = /i 2 + /i( 7 + 7*) + i(/i +p) , (3.21) 

5-i - A(3 = 2(3-i + "*7 - l*P - Pi* , (3.22) 

Ap = pp* - p{n + 7 *) - jp + 1^ , (3.23) 

Aa - 5*7 = n*a - /3* 7 - a 7 * . (3.24) 



Due to the fact that so many of the background spin coefficients are zero, the 
remaining Ricci identities give no further information. We note that we have 
made use of the vanishing of the background Weyl tensor in this calculation. 
Also for clarity we refrain from substituting for the spin coefficients from 
Eqs. ()3.2Jl -(3.7) until the very end. 
Under the (class III) rotation 

k J -> k j , V -> V , m j -> e ilp m j , m J ' -> e^rhP , (3.25) 

where A and (p are two real functions the non-zero spin coefficients transform 
as 



P -> P , 

1 .. 

7^7+ 2^ A< P ' 
e -> e + , 

(3 -> + -ie^fy . 



(3.26) 
(3.27) 

(3.28) 
(3.29) 
(3.30) 
(3.31) 
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From these transformations we see that we can in principle rotate away e — e*, 
7 — 7* and a — (3* (that is we can assume e, 7 are real and a = /?*) provided 
the function ip satisfies 



D<p = i(e-e*), A(p = 1(7 - 7*) , 6(p = i{a - /3*) . (3.32) 

However before we impose these conditions we need to check that they are 
consistent with each other and with all other equations. Suppose the first of 
these two conditions hold. Using Eq. (J3.16J) we find 

DAcp - ADcp = 2i(e~f - e* 7 *) (3.33) 

and it is straightforward to check that this agrees with one of the standard 
commutation relations. Thus we can set e — e* = and 7 — 7* = without 
imposing any extra conditions. If we assume that the third of ()3.32|) holds 
then Eq. (|3.19|) and the commutation relation involving [5, 5*} are consistent 
only if we have 

H*p* + HP = -3^0 • (3.34) 
Since p, p are real this simplifies to 

1 

2 /^P = -g/^o • (3.35) 
Replacing p here by (3.3) yields the quadratic equation 

3p 2 - 26 p + p = . (3.36) 

In order to have p real we require 

6 2 > 3p . (3.37) 

Here 9 is the expansion of the universe with 9 = 3R/R where R(t) is the scale 
factor of the general Robertson-Walker line-element and the dot indicates 
differentiation with respect to t. Using the background Einstein equations 
Gij = (po + p)uiUj + pgij, where is the Einstein tensor, we obtain the 
background Friedmann equation 

e 2 = 3p - ^ , (3.38) 

where k = 0, ±1 is the Gaussian curvature of the homogeneous space-like 
hypersurfaces t = constant. Comparing this with ()3.37j) we see that if we 
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impose (|3.34|) we cannot continue to study universes with k — +1. However 
in pQ we constructed examples of information carrying waves propagating 
through such universes and so we shall not insist on having a — f3* = 0. Thus 
from this point we assume only that we can set 7 — 7* = and e — e* = 
without loss of generality. 

We complete this section by listing (in terms of the spin coefficients) some 
derivatives of the tetrad vectors which we use extensively in the next section: 

rrii-jm 3 = (a* — (3)rrii , (3.39) 

rhi-jiji 1 = (a — P*)rrij + ((3 — a*)rhj , (3.40) 

m^jifi? = —fi*ki + pli + {(3* — a)rrii , (3-41) 

hj j = -{i + i*)h , (3.42) 

m i;j k j = , (3.43) 

///,:/' = . (3.44) 



4 The Key Result 

We now have all the equations necessary to derive our result from the wave 
equation ()2.7|) for s u using the propagation equation (|2.8J) and the divergence- 
free condition ([2.5)1 . To begin with we examine the divergence-free condition. 
Writing this condition in tetrad formalism using (j2.14|) and then contracting 
with rrii yields 

5*s = -2(a- (3*)s . (4.1) 

Thus we now have an expression for the derivative of s in the direction of 
m\ In terms of the spin coefficients we can write 

9 = ^p-3fi (4.2) 

and 

0/ = -20p . (4.3) 

Multiplying the propagation equation ()2.8j) by m l wP and using these expres- 
sions we arrive at (since tfij = —<f>ki) the following equation for the derivative 
of s in the direction of k % : 
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Da=(fi- |p) s + ±11 . (4.4) 
Similarly we multiply the wave equation (J2.7)) by m l m^ and this yields 

s^-kTUimj - h)s IJ in,in l - Qfl + -0 2 ) s + (p - s 



-IFm^- - ^011 . (4.5) 



Writing s v;fe ; & as g s^-m and using Eq. (|2.14jl we can (on substitution from 
Eqs. f)3.39j) - (|H.44j) ) write the first term of this wave equation as 



s^-kUiiUij = 2[S(5*s)-A(Ds))-2(j + j*-fM)Ds + 2s5(a-(3*) 
+ 2 09 - a*)5*s + 4(a - /T)5s + 2s<P(/3 - a*) 
- 2p*As + 8(f3-a*)(a- (3*)s-2(fi*p* + fip)s . (4.6) 

Replacing Ds here by ()4.4j) . 5*s by (j4.1j> and using Eq. (|3.19|) to write 5*a—5(3 
(and its complex conjugate) in terms of undifferentiated spin coefficients this 
equation simplifies to 



s lj ' k . k mirnj 



-2A/i + 3Ap-2(^-^ [y-~\p 

+ \w - (-0 - //) n - An . 



s + -0As 

o 



(4.7) 



We have also used Eqs. (j3.8jl . ()3.9|) and (|4.2j) to write the equation in this 
form. 

We now turn our attention to the second term of the wave equation (|4.5j) . 
As a consequence of ()2.14j) we have 



s^mirrij = s + 2smmi 



(4.i 



where the dot indicates covariant differentiation in the direction of the 4- 
velocity u*. Writing out this derivative explicitly (i.e. write s = s.^u 1 ) 
and noting that u l = l l + \k % the above equation becomes (on account of 



7-7 



e — e 



0) 



s lJ mimj = As + --Ds . 



(4.9) 
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The anisotropic stress IP- 7 satisfies the same equations as s l i (c.f. (12. 4 j) and 
()2.5|) and thus we also have 

mirrij = AU + ^DU . (4.10) 

Now reconstructing the wave equation ()4.5|) using Eqs. (|4.7)L ()4.9|) and ()4.10|) 
yields 



3Ap - 2A/i - 2 Qfl - /i) (/J - ^) - -0 - 

+ ^ + ^ ) s -^ + l^ s = -lDn-^n-p n. (4.n) 

The background Raychaudhuri equation is 

^ = -^ 2 -^(/io + 3p) . (4.12) 

Making use of this, the expressions for A/i and Ap given by Eqs. (J3.17)) 
and 1)3.21)) respectively and Eq. ()4.2)) we find that the wave equation ()4. 1 1|) 
reduces to the remarkably simple form 

(Ho + p)s= -DTl-pTl . (4.13) 

The necessity for II ^ follows from this equation because it is immedi- 
ately clear that if the perturbed anisotropic stress vanishes (<(=> II = 0) then 
the perturbed shear is also zero (•£=> s = 0) if p + P 7^ and we have no per- 
turbations (i.e. no gravitational waves). We note that the right hand side of 
equation (|4.13|) vanishes if II satisfies the differential equation UIl-l- pll = 0. 
In this case we must have s = (since po + P 7^ 0) and thus from the 
propagation equation ()4.4j) we must have II = 0. 



5 Summary 

We have investigated general properties of some gravity wave perturba- 
tions of isotropic cosmologies. Specifically we have considered gravitational 
waves carrying arbitrary information. By this we mean that the perturba- 
tions describing the waves have an arbitrary dependence on a real-valued 
function 0(x*). Our main result hinges on two important properties; in order 
to have pure gravity waves the hypersurfaces <j)(x l ) must be null and shear- 
free. Now our main result (proved in Section 4) is that under the physically 
reasonable assumption that p +P 7^ (where po and p are the proper density 
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and pressure of the cosmic fluid) the perturbations describing the gravita- 
tional waves must be accompanied by the presence of anisotropic stress in 
cosmic fluid. 
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A Notation and Basic Equations 

Some notation and tensor quantities associated with the vector field u % in 
Section 2 are required. Covariant differentiation is indicated with a semi- 
colon and covariant differentiation in the direction of u % is denoted by a dot. 
As usual square brackets denote skew-symmetrization and round brackets 
denote symmetrization. Thus the 4-acceleration of the time-like congruence 
is 
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u* := uVj u J ' . (A.l) 
With respect to u l and using the projection tensor 

hij := gij + UiUj , (A. 2) 

the covariant derivative of the 4- velocity u^j can be decomposed into 

Ui-j = Uij + Gij + -9 hij - iii uj , (A.3) 
o 

with 

Wij ■= U[i;j] +U[iUj] , (A.4) 

the vorticity tensor of the congruence, 

Oij := U{i-j) + U(iUj) -^9 , (A.5) 
the shear tensor of the congruence and 

0:=uV, (A.6) 

the expansion (or contraction) of the universe. 

The Riemann curvature tensor Rijki is defined by the Ricci identities 

Ui-ik - Ui. k i = Rijkiu J , (A. 7) 

and Einstein's field equations take the form 

Rij ~ -j9ij R — Tij , (A. 8) 

where Rij := R k ikj are the components of the Ricci tensor, R := R\ is the 
Ricci scalar and is the symmetric energy-momentum-stress tensor. We 
note that the coupling constant has been absorbed into T^. With respect 
to the 4-velocity field u % the energy-momentum-stress tensor can be decom- 
posed as 

T ij = fiQ u l u j + p h ij + q i u j + q j %j + n ij , (A. 9) 

with 

q i u i = , ir ij Uj = , 7T*j — , (A.10) 

and = 7T-' 1 . Then /xq is interpreted as the energy density measured by 
an observer with 4- velocity u\ q l is the energy flow (such as heat flow) 
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measured by this observer, p is the isotropic pressure and is the trace- 
free anisotropic stress (due for example to viscosity). 

With respect to u l the Weyl tensor may be decomposed into its "electric" 
and "magnetic" components given respectively by 

Eij = C ikj i u k u l , H i:j = *C ikjl u k u l , ( A. 1 1 ) 

where *C^j\ = | r] ik mn C mn ji is the dual of the Weyl tensor (the left and 
right duals being equal), r] m = y/=g e ijkl with g = det(^) and e ijkl is the 
Levi-Civita permutation symbol. The explicit expression for the Weyl tensor 
in terms of Ey and Hy is not needed here but can be found in 
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